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ABSTRACT 



Context. The colliding-wind binary 77 Carinae exhibits soft X-ray thermal emission that varies strongly around the periastron passage. 
It has been found to have non-thermal emission, thanks to its detection in hard X-rays using INTEGRAL and Suzaku, and also in y-rays 
with AGILE and Fermi. 

Aims. This paper attempts to definitively identify r\ Carinae as the source of the hard X-ray emission, to examine how changes in the 
2-10 keV band influence changes in the hard X-ray band, and to understand more clearly the mechanisms producing the non- thermal 
emission using new INTEGRAL observations obtained close to periastron passage. 

Methods. To strengthen the identification of rj Carinae with the hard X-ray source, a long Chandra observation encompassing the 
INTEGRAL/1SGR1 error circle was analysed, and all other soft X-ray sources (including the outer shell of r\ Carinae itself) were 
discarded as likely counter-parts. To expand the knowledge of the physical processes governing the X-ray lightcurve, new hard X-ray 
images of r\ Carinae were studied close to periastron, and compared to previous observations far from periastron. 
Results. The INTEGRAL component, when represented by a power law (with a photon index T of 1.8), would produce more emission 
in the Chandra band than observed from any point source in the ISGRI error circle apart from r\ Carinae, as long as the hydrogen 
column density to the ISGRI source is lower than A% g 10 24 cm 4 . Sources with such a high absorption are very rare, thus the hard X- 
ray emission is very likely to be associated with r\ Carinae. The eventual contribution of the outer shell to the non-thermal component 
also remains fairly limited. Close to periastron passage, a 3-cr detection is achieved for the hard X-ray emission of 77 Carinae, with a 
flux similar to the average value far from periastron. 

Conclusions. Assuming a single absorption component for both the thermal and non-thermal sources, this 3-cr detection can be 
explained with a hydrogen column density that does not exceed N H < 6 x 10 23 cirr 2 without resorting to an intrinsic increase in the 
hard X-ray emission. The energy injected in hard X-rays (averaged over a month timescale) appears to be rather constant at least as 
close as a few stellar radii, well within the acceleration region of the wind. 
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1. Introduction 



Located in the Carina nebula at a distance of 2.3 kpc (Smith 
2006 ), 77 Carinae is not only one of the most massive objects of 
our Galaxy (80-1 20 M ; |Davidson & Humphrey s| 1997] |Hillier 
et al.||206T| l, but also one of the most peculiar and impressive. 
During its "Great Eruption", it briefly became the second bright- 
est object in the sky in 1843, before fading to a V magnitude my 
of 8 . A second notable eruption took place in 1 890. In the middle 
of the XXth century, it started to slowly and irregularly brighten 
again , up to its current value of my ^ 5 (see e.g. Davidson et al. 
T^99||Fernandez-Lajus et al.|2Q10] >. 



X-rays: binaries - X-rays: individuals: rj Car - X-rays: individuals: IE 1048.1-5937 - X-rays: 



the object makes observations in the optical and ultra-violet en- 
ergy ranges very difficult, rj Car continues to eject matter through 
energetic stella r winds, with an estimated mass-loss rate of 10~ 4 - 
10 3 M Q yr 4 (jAndriesse et al.|l978l|Hilrier et al.|200lj |Pittard 



The extended bipolar nebula that can be seen around 77 Car, 
the so-called homunculus, is due to the huge quantities of matter 
(10-20 M ; |Smifh et al.|2Q03| > ejected during the 1840s, while 
the little homunculus is due to the ejection of approximately 
0.1 Mq during the second burst in the 1890s (Ishibashi et al. 
2003 Smith 2005). The circumstellar gas that now enshrouds 
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|& Corcoran|2002[|van Boekel et al.|2003| l. 



1.1. t) Carinae as a binary system 

Although many questions remain open, there is now strong ev- 
idence that rj Car is a binary system. For example, a period 
of ^ 5.5 yr, which has been stable over decades, was inferred 
from observations in the radio ( |Duncan et al.|1993 1, millimeter 
(mm; |Abraham et al.||2005| l, optic al (|Damineli||1996l |Damineli| 
l et al.|2000|>, nea r-infrared (n ear-IR;|Whitelock et al.|1994[|2004t 
Damineli 1996 ), and X-ray ( |Corcoran|2005| l domains. |Damineli| 
|et al.| ( |2008b| l found an orbital period of ~ 2022.7 + 1.3 days 
(5.53 yr) by combining radial velocity variations from multiple 
lines in optical spectra along with broad-band IR, optical, and 
X-ray observations. 

The primary star is most likely a luminous blue variable 
(LBV; Davids on & Humphreys|1997] ), a short phase during the 
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evolution of massive stars that occurs after leaving the main se- 
quence and leads to the object becoming a Wolf-Rayet star. The 
secondary star is probably a late-type nitrogen-rich O or WR 
type ( jlping et al.|2005||Verner et al.|2005] ). |Mehner etaL] < [20T0| > 
derived a zone in the Hertzsprung-Russel diagram where the sec- 
ondary is most likely to lie (see their Fig. 12). The orbit has a 
semi-major axis of 16.64 AU (jHillier et al. 2001) and is very 
eccentric (e =s 0.9; Nie lsen et al.|2007) >. These values lead to a 
periastron distance r per i astron of 1 .66 AU, which should be com- 
pared to the (poorly known) primary star's radius R+,u believed 
to be ^ 0.7 AU ( |Corcoran & H amaguchi 2007), but which might 
be as high as 1 AU (Damine ti|1996| ). 

The orientation of the binary system remains unclear. A 
number of studies have concluded that the secondary moves be- 
hind at periastron (|Damineii|1996"}|Pittard et al.|1998[|Corcoran 
et al.|2001 ||Corcoran|2005||Akashi et al.|2006||Hamaguchi et al. 
20071 INielsen et al. 112007 1 Kashi & Soker|2008; Okazaki et al. 



2008 ; Par kin et al.|2009[), but others seem to favour the opposite 



situation (jFalceta-Gongalves et al.|2005j|Kashi & Soker 



2007) 



Smit 



while quadrature has also been suggested (Ishibashi 2001 
|et al.|2004) . 

1.2. X-ray observations 

Observations of 77 Car in the X-ray domain are relatively unaf- 
fected by absorption compared to most other wavelengths, thus 



allowing an investigation deep into the system (see Corcoran & 
Hamaguchi|2007l ). 

First shown by the Einstein satellite (Chlebowski et al. 



1984), and later confirmed by Chandra observations (Seward 



et al.|200T 1, the X-ray emission of 77 Car has two distinct com- 
ponents: 

- a softer, spatially extended (up to about 20"), and inhomo- 
geneous thermal (kTsx — 0.5 keV) X-ray component (named 
77 SX), which dominates the spectrum mostly below 1.5 keV, and 
is probably associated with the interaction between the stellar 
wind and the interstellar matter; 

- a harder, point-like thermal (kTux — 4.7 keV) X-ray compo- 
nent (77 HX), centered on the stellar system, that dominates the 
spectrum in the 2-10 keV range and is probably linked to the 
wind collision between the two massive stars that form the bi- 
nary system. 

The X-ray spectrum is coherent with a colliding- wind binary 
(CWB) interpretation (|Usov|1992||Stevens et al.|1 992 1 Pitt ard&| 
Stevens|19971|Pittard et al.|19981|Ishibashi et al.|l999l|Pittard &| 
Corcoran|2002t|Corcoran|2005||Pittard|2007||Parkin etal. 12009^ 7 



In this framework, the slow and dense stellar wind from the LBV 
500 km/s, Mi - 2.5 x 10~ 4 M /yr; 



(v, 



00,1 
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2002 ) collides with the higher- velocity, lower-density wind from 
the hotter and luminous companion star (vqq ? ^ 3000 km/s, 
M 2 - 1.0 X 10~ 5 M /yr; |Pittard & Corcoran |2002| . The X-ray 
emitting region where the hydrodynamical shock takes place is 
called the wind-wind collision region (WCR). 

The X-ray lightcurve of 77 Carinae is of particular interest 
: it exhibits impressive variations, with the orbital periodicity, 
and undergoes an extremely deep and relatively long minimum 
close to periastron (see Sect. [2] for details). Variations at many 
other wavelengths are also coincident with the periastron pas- 
sage; for example, optical excitation lines undergo a rapid de- 
crease in intensity (Damineli et al.|2008a|l. Although the X-ray 



emission represents only a tiny fraction of the bolometric lumi- 
nosity (Lx/^boi — 10~ 7 ), understanding the origin of the hard X- 
ray emission is crucial, as it is compact and not affected by cir- 
cumstellar absorption. Moreover, the X-ray lightcurve exhibits 



short-term flares, with a quasi-period of 84 d (jlshibashi et al.| 
|T999l|Corcoran|2005| l. poffat & Corco"ran| ( |2009| ) discuss differ- 
ent interpretations of these flares, favouring an explanation based 
on the presence of clumps in the wind of the LBV. 

A hard X-ray tail (i.e., at a flux higher than an extrapola- 
tion of the thermal X-ray component 77 HX) was once detected 
in the 77 Car region by the PDS instrument onboard BeppoSAX 
( Viotti et al.|2004|l. Since then, it wa s unambiguously confirmed 
by INTEGRAL ( Leyder et al.|2008) , and reobserved by Suzaku 
( Sekiguchi et al.|2009| l. This hard X-ray detection is important, 
as it proves that non-thermal particle acceleration operates in the 
wind collision, and implies that gamma-ray emission is likely in 
the MeV and GeV energy ranges. The emission mechanism is 
believed to be inverse Compton (IC) scattering of UV photons 
by electrons accelerated in the hydrodynamical shock between 
the stellar winds ( Leyder et al.|2008" 1. 



1.3. y-ray observations 

In the y-ray domain, the Carina region has been observed since 
July 2007, in i tially by the AGILE satellite dTavani et al.|2008) . 
Tavani et al.] p009 ) report the detection of 1AGL J1043-5931 



above 100 MeV, based on observations taken with the Gamma- 
Ray Imaging Det ector instrument (GRID; 30 MeV-30 GeV; 
|Prest et aI7||2003) l. The flux of this source, averaged over all 
the observations, is (3.7 ± 0.5) X 10 -7 photon cirT 2 s~' (above 
100 MeV). |Tavani et al.| ( |2009| l claim that this source, located at 
(/, b) = (287.6, -0.7) (with an error of 0.4°), is probably asso- 
ciated with 77 Car (located 0.07° away). They also report that a 
2-day y-ray flare was observed in October 2008. 

In the same high-energy range, the Carina region has also 
been monitored since August 2008 by the Large Area Telescope 
(LAT; |Atwood et al.|2009] l instrument onboard the Fermi satel- 
lite. According to the Fermi LAT first source catalog ( Abdo et al. 
|20T0| ), the source 1FGL J 1045 .2-5942 is located at a position 
of (l,b) = (287.6252,-0.6388) (with a 95%-confidence semi- 
major axis of 1.4'). This position is 1.7' away from 77 Car, and 
the sources cannot be formally associated with each other. The 
flux is (2.6 + 0.2)x 10- 10 erg cm" 2 s _1 (in the 100 MeV-100 GeV 



energy range). In both this catalog and Tavani et al. (2009 1, it is 
indicated that there might be an association between the AGILE 
and the Fermi sources (although they are separated by 4'). 

In this paper, new hard X-ray observations performed with 
INTEGRAL close to the periastron passage of 77 Carinae are pre- 
sented. Section |2"1de scribes the behaviour of the X-ray lightcurve 
of 77 Car. Section|3|deals with the data reduction and analysis of 
the new INTEGRAL observations; while Sect. |4] describes the 
newly obtained images and spectra with respect to the previ- 
ous observations in the hard X-ray domain. Section[5]is focused 
on definitively identifying the counter-part of the INTEGRAL 
source as 77 Car. Section|6]discusses these new observations, and 
conclusions are finally drawn in Sect. [7] 

2. The X-ray lightcurve 

2.1. RXTE observations 

Observations obtained by the Rossi X-ray Timing Explorer 
(RXTE; Swank 2006 1 since 1996 cover over two orbital cycles 
of 77 Carinae, including the three minima from 1998, 2003.5, and 
2009. The X-ray lightcurve around the last X-ray minimum is 
shown in Fig . [T|(see |Corcoran|2005 | and Moffat & Corcoranf2009 
for details about RXTE data analysis), with superimposed red 
lines indicating the new INTEGRAL observations around perias- 
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tron. The blue triangles (for which the axis is given on the right 
side) are the absorption column density values obtained from 
Chandra and XMM-Newton observations around the 2003 mini- 
mum, as reported by Hamaguc hi et al.| ([2007 ) in their Table 4. 
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phase ($) 
3 



3.05 




54700 
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Fig. 1. RXTE lightcurve of 77 Carinae (in the 2-10 keV energy 
range), from August 2008 to June 2009, corrected for instru- 
mental background (from Corcoran et al. 2010, in prep.). The 
bottom axis gives the time in MJD, while the top axis gives the 
phase of rj Car (based on the ephemeris from Corcoran 2005 



The start of the minimum is close to the periastron passage from 
2009 (O = 3). Red lines indicate the new INTEGRAL observa- 
tions, obtained just after the X-ray minimum. Blue triangles are 
the Nh measured during the previous minimum (with their axis 
on the right side), as reported by Hamaguchi et al. { 2001) . 



The general behaviour of the X-ray lightcurve is closely de- 



scribed by the CWB model (Pittard et al. 1998 


Ishibashi et al. 


|1999||Pittard & Corcoran||2002| |Corcoran||2005 


. Although the 



X-ray minimum {i.e. phase O = 0) probably occurs around the 
time of the periastron passage, both events might be separated 
by up to a few months (Da mineli et al.||2000[ |Corcoran et al.| 
|2001| l. The scenario is as follows: when the two stars are far 
from each other {i.e. 3> > 0.1), the X-ray flux is roughly con- 
stant. After apastron (cD > 0.5), the secondary approaches the 
primary, the wind interactions become more important, and the 
X-ray flux increases steadily over approximately 2 years. Then, 
close to periastron, the X-ray flux drops sharply by a factor of 
more than ^ 100 (as is shown by an XMM-Newton observation 
taken during the minimum; Hamaguchi et al. 2007} for a dura- 



tion of 80-90 days. The key physical processes explaining this 
sudden decrease remain uncertain. After remaining at this mini- 
mum level for approximately 3 months, the X-ray flux returns to 
its roughly constant value as the stars recede from each other, be- 
fore starting a new cycle. The X-ray lightcurve follows the same 
periodicity of 2024 + 2 days as observed at other wavelengths 
Pshibashi et al.|1999||Corcoran|2005||Hamaguchi et al.|2007| ). 



2.2. Interpretation of the X-ray minimum 

To explain the minimum in the X-ray emission close to peri- 
astron, different scenarios have been proposed (not necessarily 
mutually exclusive). There are two broad categories of proposed 
models: some are based on a change in the intrinsic emission 
(e.g. a modification of the WCR itself), while others postulate 
a variation in the observed emission only (either by occultation, 
or by any process leading to an increase in the column density). 
These two types are developed below. 

The first category is based on the fading (or complete disap- 
pearance) of the WCR. In this model, as the separation between 
the stars decreases near periastron, they become so close that the 
WCR enters into the acceleration region of the secondary's stel- 
lar wind. The momentum balance becomes unstable, the primary 
star's wind overwhelms the other one, and the WCR collapses 
onto the surface of the secondary star. 

The second explanation of the decrease in the observed emis- 
sion is based on the eclipse model, in which the WCR becomes 
totally occulted around periastron. This can be caused by either 
the primary star itself, or its optically-thick wind; both cases re- 
quire the secondary to move behind the primary. However, the 
minimum is too long to be explained only by an eclipse due 
solely to the star itself, and would not be constant for such a 
long period because the stellar motion close to periastron is very 
rapid ([Hamaguc hi et al.|2007| >. 

Another explanation of the observed X-ray variation is a 
temporary increase in the mass-loss rate. In this hypothesis, 
an increase in the mass-loss rate leads to more material be- 
ing ejected into the interstellar medium (ISM), and therefore a 
shielding of the UV photons, thus encouraging the formation of 
dust and an increase in absorption. 

The X-ray minimum can therefore be explained by either 
a decrease in the plasma emission measure (EM), an increase 
in the column density, or a combination of those two options. 
The column density A^h - 5 x 10 24 cirT 2 necessary to reduce 
the X-ray flux by two orders of magnitude from the maximum 
to the minimum is ind eed much higher than th e value observed 
{Nh - 5 x 10 23 cm" 2 ; Hamaguchi et al.||2007 1, hence the need 



for an additional explanation, distinct from the column density 
variation. 

Hard X-ray observations at the X-ray minimum, such as 
those provided by INTEGRAL, which are mostly unaffected by 
absorption A^h 5 10 24 cm" 2 , can help distinguishing between 
these different models. 



3. INTEGRAL observations and data analysis 

The ESA INTEGRAL y-ray satellite carries 3 coaligned instru- 
ments dedicated to the observation of the high-energy sky, from 
3 keV up to 10 MeV, along with an optical monitoring camera 
( |Winkler et al.||2003[). The INTEGRAL soft gamma-ray imager 
(ISGRI; Lebrun et al. 2003 ), the most sensitive imaging detec- 
tor between 15 and 200 keV currently flying, offered the first 
apparent detection of r\ Car at hard X-rays ( Leyder et al.|2 008 ) 



The new INTEGRAL observations are indicated by the red 
lines superimposed on the RXTE lightcurve in Fig.[T] They were 
scheduled for a period covering 80 days after the expected 2009 
minimum, as the previous minima from 1998 and 2003.5 lasted 
80-90 days. Unfortunately, the 2009 minimum happened to be 
half as long as the previous two X-ray minima, explaining why 
only a fraction of the INTEGRAL observations correspond to the 
minimum. 
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Two images were produced, to compare the behaviour of 
77 Car close to and away from periastron. The first image contains 
all available data located within 15° of 77 Car, and which covers 
the minima; this includes data from the 2003.5 X-ray minimum, 
and from the 2009 X-ray minimum (see Table[T]l. The second im- 
age was compiled using all available public data located within 
15° of 77 Car and not included in the first image; this second 
image is thus comparable to the one from Leyder et al. (2008), 
although the effective exposure time is significantly longer. A to- 
tal of 362 (respectively 2027) pointings were added to produce 
the first (respectively second) image, for a deadtime-corrected 
good exposure of 1.1 Ms (respectively 5.5 Ms). ISGRI point- 
ing sky images were produced with a pre-release of OS A[J ver- 
sion 9, using standard parameters. The image-cleaning step used 
an input catalogue of all INTEGRAL sources with fixed source 
positions, and was applied independently of the source strength 
(thus allowing for negative source models) to avoid introducing 
any biases during the process. 



Table 1. INTEGRAL data available for 77 Car close to the peri- 
astron passage; only data coincident with the X-ray minimum 
were selected {i.e. phase < <t> < 0.04). 



Minimum 


Revolutions 


Time [MJD] 


Phase O of 77 Car 


2003.5 


88-91 


52 824-52 833 


2.000-2.005 


2009 


770-788 


54 865-54920 


3.009-3.036 



Broad-band (17-80 keV) mosaic images of the field were 
produced, along with narrower bands in 12 energy ranges. The 
final mosaic images were compiled using equatorial coordinates 
with a tangential projection, using an oversampling factor of two 
with respect to the individual sky images. The photometric in- 
tegrity and accurate astrometry were ensured by calculating the 
intersection between input and output pixels, and weighting the 
count rates with the overlapping area. 



4. INTEGRAL images and spectra 



The previous analysis of ISGRI data outside periastron (Leyder 
|et al.|2008] l identified three sources in the field of 77 Car: 77 Car 
itself, IGR J 10447-6027, and the anomalous X-ray pulsar (AXP) 
IE 1048.1-5937. 

The new ISGRI 17-80 keV image close to periastron only 
(i.e., the first image, shown in the left panel of Fig. [2} indicates 
that the hard X-ray emission of 77 Car close to the X-ray mini- 
mum is around 0.16 ± 0.05 cnts -1 . Although neither 77 Car nor 
IGR J10447-6027 are detected at a level of 5-cr, excesses are 
clearly observed at both positions at a level of 3-cr. 

The new image far from periastron (i.e. the second image, 
shown in the right panel of Fig. [2]) indicates that both 77 Car 
and IGR J 10447 -6027 are very well detected, providing more 



accurate positions and average fluxes than Leyder et al. (2008). 
IGR J 10447 -6027 is thus likely to be a persistent source. In con- 
trast, the AXP IE 1048.1-5937 is no longer clearly detected (al- 
though an insignificant excess is observed at its position), as op- 
posed to the situation reported in |Leyder~e t al. (2008), and is 
therefore probably variable. 

Table [2] summarizes the best-fit position, error circle radius, 
intensity (or 5-cr upper limit), and significance for the sources 




Fig. 2. ISGRI significance image of 77 Car close to periastron 
(on the left) and outside periastron (on the right). Both images 
are in equatorial coordinates, cover the same energy range (17- 
80 keV), and use the same colourmap (significance goes from 
2<x to 7<x). The positions of 77 Car and IE 1048.1-5937 are indi- 
cated by the white and black crosses. The white and black circles 
represent the best-fit positions listed in Table [2] 

detected outside periastron; it also indicates the intensity and 
significance at periastron. 

The new ISGRI best-fit position, now only 1 .6' away, is sig- 
nificantly closer to 77 Car than the previous one ( |Leyder et al. 
2008), and the error circle is slightly smaller, decreasing to 2.4'. 
These new values strengthen the association of the INTEGRAL 
source with 7/ Car. 

The average source flux for 77 Car observed by ISGRI (ex- 
tracted from the second image by fixing the PSF to 6' and 
the source to its optical position) is 0.18 ± 0.02 count s in 
the 17-80 keV energy range. This corresponds to a flux of 
1.2 x 10~ 11 ergcm _2 s _1 when assuming a photon index of 1.8 
(which leads to the best fit of the data, as will be shown below). 



Based on a distance of 2.3 kpc (Smith 2006), the hard X-ray 
luminosity is 8 x 10 33 ergs -1 (17-80 keV). 

The limited number of points in the ISGRI spectrum (ex- 
tracted from data obtained far from periastron) only allows sim- 
ple models to be fitted to the hard X-ray spectrum. A power 
law alone gives a photon index T of 1.7 ± 0.5. As there are no 
soft X-ray data simultaneous with the INTEGRAL observation, 
an archival BeppoSAX observation (from December 29, 1996) 
was used with the ISGRI data. Note however that this X-ray 
observation is adopted as a reference only, and might slightly 
overestimate the contribution of the thermal component in the 
ISGRI band. The model applied was a combination of an opti- 
cally thin, thermally emitting plasma (Me we et al.| [T985 1986; 
Kaas trap992 Liedahl et al.||1995] ) with a power law, both un- 
dergoing the same amount of interstellar photoelectric absorp- 
tion (using cross-sections from |Morrison & McC ammon 1983 ) 
and Thomson scattering (i.e. wabs*cabs*(mekal+powerlaw) 
in Xspec terms). The parameters derived are a photon index of 
F = 1.8 + 0.4 and a temperature of kT = 4.3 keV (giving a re- 
duced^- 2 of 3.6 for 78 degrees of freedom, down to 1.5 when 
assuming a 5% systematic error). 

This new value for the photon index is in closer agree- 
ment with Suzaku observations than the previous value based 
on INTEGRAL data (|Leyder et alj|2"0"08"l). The contribution of 



1 The Offline Scientific Analysis (OSA) software is available from 
the ISDC website : http : // www . i sdc . uni ge . ch 



the optically thin thermal plasma emission (mekal) in the first 
ISGRI bin remains important (about 40%), and is poorly known 
because of the lack of a simultaneous soft X-ray observation 
(Hamaguchi et al. 2007 report variations by a factor 2 above 
3 keV). 
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Table 2. Sources detected around 77 Car in ISGRI mosaics. 



Source 



Intensity [cnts ] 



Position (J2000) 



Distance ['] Error circle radius ['] Significance 



Close to periastron (Fig. 2 left panel) 



77 Car 
IGR J10447-6027 



0.16: 
0.20: 



0.05 
0.05 



Far from periastron (Fig. PI right panel) 

. . — ~ W , : 



3.3 
4.2 



RA = 10 E 44 m 57 s , Dec = -59°42'27" 
RA = 10 h 50 m 14\ Dec = -59°53'10" 
RA = 10 h 44 m 36'\ Dec = -60°25'50" 



77 Car 
IE 1048.1-5937 
IGR J10447-6027 



0.18 ±0.02 

< 0.11 
0.14 ±0.02 



1.6 
0.6 
1.9 



2.4 
6.2 
3.0 



8.4 
2.9 
6.3 



Notes. The flux and significance values were extracted from the 17-80 keV energy band with the best-fit position, and with the PSF size set to 6'; 
the values close to periastron are for the 3-cr detections; the upper-limit is at the 5-cr level; the distance refers to the extent between the positions 
derived from the ISGRI observations and the published positions; the error circle corresponds to a 90% probability. To extract the values from the 
observations close to periastron, the position of the sources were fixed to their best-fit positions derived from the observations far from periastron 
(since the lower significance of the former implies that they provide far less precise positions). 



The quality of the fit can be improved by allowing the chem- 
ical abundance of the plasma to vary and by ignoring the softest 
part of the BeppoSAX spectrum (which would require the addi- 
tion of another optically thin thermal plasma model, i.e. mekal, 
emitting at a different temperature, to be correctl y modelled), 



Viotti et al. 



leading to a reduced x 1 of about 2 (for details, see 
2004); this however does not affect the values derived for the 
photon index. 

The unfolded spectrum outside periastron is shown in Fig. [3] 
The MECS observation from December 29, 1996 (at phase <t> = 
0.83; black points) is only meant to provide a reference spectrum 
of the source outside periastron. The red points are the ISGRI 
spectrum far from periastron. The green point is the ISGRI 3-cr 
detection close to periastron. The unfolded model is shown in 
magenta, the thermal component being plotted in dark blue and 
the power law in light blue. The fluxes observed close to and 
outside periastron are consistent. 
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Fig. 3. BeppoSAX/MECS (black) and INTEGRAL/lSGRl (red) 
unfolded spectra of 77 Carinae outside periastron (the last ISGRI 
point is a 5-cr upper-limit); ISGRI 3-cr detection close to peri- 
astron (green). The unfolded model is the magenta long-dashed 
line; with the thermal and non-thermal components respectively 
shown by the dotted dark blue and light blue lines. 



Finally, the data were also fitted with a model 
in which both thermal and non-thermal compo- 
nents are subject to different absorption columns 
(wabs*cabs*mekal+wabs*cabs*poweiTaw). The best-fit 



solution was obtained with values very similar to those found 
when they are tied together, although the absorption on the 
power law is poorly constrained. However, hydrogen column 



densities higher than Nh £1.5 xlO cm can be clearly ruled 
out, as the hard X-rays become so absorbed that the first ISGRI 
bin can no longer be correctly modelled, independently of the 
model chosen. 



5. Counterpart of the non-thermal source 



Leyder et aL] (|2008) associated the INTEGRAL hard X-ray de- 
tection with 77 Carinae, based mostly on this source being by far 
the most X-ray luminous in the immediate vicinity. However, all 
other X-ray sources within the ISGRI error box (based on a long 
Chandra observation) are very unlikely to be the counterpart of 
the ISGRI source, as shown below. 

To compare the hard X-ray spectrum to all X-ray sources 
located within the ISGRI error circle, the longest imaging obser- 
vation of the 77 Car field taken by Chandra was chosen. It is an 
88.2 ks observation taken in August 2006 (ObsID is 6402). This 
observation was taken far from periastron, which is desirable. 
The data were reduced using CIAC0 version 4.1.2, with stan- 
dard parameters and methods. After data reduction, application 
of the appropriate good time intervals (GTIs), and selection of 
the energy range of interest (0.3-10.0 keV), the exposure time 
was reduced to 84.9 ks. 

From a Chandra image in the 5-10 keV energy range, all 
sources with more than a few counts and located inside or close 
to the ISGRI emission were selected (see Fig. |4j the ISGRI er- 
ror circle is drawn in red, while the concentric 4.5' selection re- 
gion is shown in blue; sources are labelled in order of decreasing 
brightness, provided that they have more than 7 counts). Events 
linked to either 77 Car or its nebula were not considered, i.e. those 
within 1.5' of the source position or in the black box (which is 
an artifact due to the huge brightness of 77 Car); both exclusion 
regions are indicated in black in Fig.|4] In this image, 77 Car is by 
far the brightest object, with significantly more than ten thousand 
counts inside the black circle. The following sources (numbered 
1, 2, and 3) contain, respectively, about 100, 50, and 40 counts. 

It is known that the brightest sources in this field are associ- 
ated with massive stars, while the majority of fainter sources are 
believed to be pre-main sequence (PMS) stars (e.g. Antokhin 
et al.||2Q08] >. Those sources are most accurately represented by 



2 The Chandra Interactive Analysis of Observations (CIAO) soft- 
ware is available from the Chandra X-ray Center (CXC) website : 
http : // cxc . harvard . edu/ciao/ 
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Fig. 4. X-ray sources detected in the Chandra 5-10 keV image 
(limits range from 3 to 10 counts). The ISGRI error circle is 
shown in red for reference, the C/ianc/ra-selected sources are lo- 
cated up to 4.5' from the ISGRI position (within the blue circle), 
and the exclusion regions (a 1.5' circle around 77 Carinae plus a 
box to avoid including an artifact) are drawn in black. 
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Fig. 5. Chandra photon counts expected in the 5-10 keV energy 
range as a function of the hydrogen column density (Nh), for a 
photon index of T=1.4 (bottom curve), 1.8 (middle curve), and 
2.2 (top curve). These extrapolations of the ISGRI observations 
are based on a power-law model that includes interstellar photo- 
electric absorption (wabs-powerlaw). 



an optically thin thermal plasma emission model taking into 
account the interstellar photoelectric absorption (wabs*mekal), 
thus they have hardly any contribution to the hard X-ray range. 

To estimate the counts expected in the Chandra image from 
the ISGRI source, its spectrum was fitted with a simple absorbed 
power-law model including interstellar photoelectric absorption 
(wabs*powerlaw) and extrapolated down to the Chandra en- 
ergy range (see Sect. [4j. The Chandra source counts expected 
in the 5-10 keV energy range are shown in Fig. [5] The three 
curves correspond to three different values of the photon index 
(T=1.4, 1.8, and 2.2), thus covering the range of possible values 
as obtained from the ISGRI observations. 

Adopting the lower limit from the range of hydrogen col- 
umn densities found by Hamaguchi et al. (2007] Nh 



40x1 22 cm -2 ), because this lower value is observed during most 
of the orbit, the number of counts expected in the 5-10 keV en- 
ergy range from the extrapolation of the INTEGRAL spectrum is 
of the order of 4300. This is thus nearly a factor of three smaller 
than the number of counts originating from 77 Car detected with 
Chandra in the same energy range. In addition, as the 2-10 keV 
flux in this Chandra observation is about 5 times lower than the 
one from the BeppoSAX observation shown in Fig. [3] the ratio 
of the thermal model fitting the MECS data to the power-law 
extrapolation from the ISGRI data is understood. 

No source other than 77 Car detected with Chandra within 
the ISGRI error circle has more than approximately 100 counts 
in the 5-10 keV energy range. According to Fig. [5] this means 
that if one of these other sources were the counterpart to the 
ISGRI emission, the hydrogen column density towards that 
source would have to be at least as high as Nh <: 10 24 cm" 2 
to be compatible with both the INTEGRAL and Chandra data. 
This is a lower limit, as a source with a soft excess or a leaky 



absorber would require an even higher hydrogen column density 
to match the observations. 

Such hard X-ray emission affected by a large absorption 
could be indicative of either a Compton-thick AGN, or a highly 
obscured high-mass X-ray binary (HMXB). For the latter, out 
of the more than 250 sources lying within the Galactic plane 
(see the fourth ISGRI catalog from |Bird et al. |2010| >, only one 



HMXB has an absorption of Nh >, 10 z * cm" / . As the ISGRI 
spectrum extracted at the location of 77 Car does not match the 
typical spectrum of an accreting pulsar (i.e. a power law with a 
cutoff at ~10 keV), this is an unlikely possibility. 

It is also unlikely that a Compton-thick AGN might be the 
counterpart to the hard X-ray detection. Paltani et al. (2008) 
studied the distribution of AGNs around the quasar 3C 273, for 
which the exposure time is comparable to that available on 77 Car. 
From that sample, they derived a logA^-logS distribution for 
AGNs with an absorption column density Nh > 10 22 cirT 2 (see 
the blue curve in the bottom panel of their Fig. 9). They also 
placed an upper-limit of 24% on the number of Compton-thick 
AGNs present in their sample. Given that the 20-60 keV flux for 
77 Car is around 8 x 10~ 12 erg crrT 2 s _I (when assuming a photon 
index of 1.8) and that the size of the error circle on 77 Car is 2.4', 
the probability of having a Compton-thick AGN coincident by 
chance with 77 Car is strictly lower than 0.002%. 

An alternative explanation to the IC emission from the CWB 
was suggested by Ohm et aT] ( |2010| l. In their scenario, the ac- 
celeration of the non-thermal particles is due to the outer shell 
surrounding 77 Car, which is detectable mostly in the 0.1-1 keV 
energy range (Sew ard et al.|200T) . However, based on the long 
Chandra observation, the number of counts detected in an annu- 
lus covering this shell (and excluding a point source) is of the 
order of 750 counts (in the 5-10 keV energy range). This value 
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is an upper limit to the shell emission, as there remains some 
contribution (estimated to be at least 10%) from the wings of 
the PSF of 77 Car itself in that extraction region. Although there 
is indeed some small contribution from this shell to the X-ray 
spectrum, this is at least a factor of six lower than the number of 
counts expected from the non-thermal extrapolation from ISGRI 
data down to the same energy domain. Thus, the contribution of 
the shell to the hard X-ray emission, if any, remains fairly lim- 
ited. 

In summary, the most likely counterpart to the observed hard 
X-ray emission is 77 Car, as other sources detected with Chandra 
within the ISGRI error circle seem unable or very unlikely to 
match the flux level detected with INTEGRAL. 



6. Discussion 

Thanks to its excellent spatial resolution, INTEGRAL has pro- 
vided the strongest evidence that the non-thermal emission 
comes from 77 Carinae. All sources detected by Chandra close to 
the ISGRI emission can be ruled out as likely counterparts, and 
the probability of a Compton-thick AGN being serendipitously 
located at the same position as Eta Car is less than 0.002%. 

The general form of the X-ray lightcurve of 77 Car close to 
periastron can be explained by an increase in hydrogen column 
density, a decline in intrinsic emission, or a combination of both. 
The motivation behind acquiring these new INTEGRAL obser- 
vations was to place constraints on the relative contributions of 
these two mechanisms. 

The fraction of particles accelerated to high energies and in- 
jected into the hard X-ray tail is assumed to be constant or lower 
around periastron than during the remainder of the orbit. This is a 
reasonable hypothesis, since although the shock will be stronger 
closer to the star, the Compton cooling will be far more impor- 
tant because of the higher density of UV photons. This means 
that the electrons are cooled more rapidly, leading to a softer dis- 
tribution and a reduced emissivity. Therefore, the expected high- 
energy emission is lower near periastron. Moreover, close to pe- 
riastron, the stars are so close to each other with respect to the 
stellar radii that the winds may not have time to reach their full 
terminal speed, thus mitigating the effect of the stronger shock. 

The ISGRI 3-cr detection close to periastron shows a flux 
level consistent with the observations far from periastron. This 
means that the hard X-ray flux may not be lower near periastron 
than outside, and might even have to be intrinsically higher if the 
hydrogen column density effectively increased strongly during 
the INTEGRAL observations. 

During the spectral modelling, it was assumed that both ther- 
mal and non-thermal components are affected by the same ab- 
sorption. |LeyderetaL]( 2008 ) suggested that the emission mech- 
anism explaining the hard X-ray tail is IC emission of UV pho- 
tons by electrons accelerated in the hydrodynamical shock that 
forms between the stellar winds. If this is indeed the case, the 
hard X-rays are then formed in the same region as the soft X- 
rays, thus it is reasonable to assume that the hydrogen column 
density will be the same for both thermal and non-thermal com- 
ponents of the model. Moreover, the WCR is spatially extended 
and fairly non-uniform, thus making it difficult to mask the hard 
X-ray emission while allowing a much smaller absorption for 
the soft X-rays. 

Simulations to estimate the 17-80 keV flux were per- 
formed for a large range of the hydrogen column density 
(Nh - 1-1000 xlO 22 cm" 2 ). The model used was a com- 
bination of an optically thin thermal plasma emission with 



a power law, both affected by the same amount of in- 
terstellar photoelectric absorption and Thompson scattering 
(wabs*cabs* (mekal+powerlaw)). The temperature was fixed 
to kT = 5.1 keV and the photon index to F = 1.8. In addition, the 
model was scaled to match either the average 2-10 keV flux ob- 
served by RXTE far from periastron 7 x 10~ n erg crrT 2 s _1 ), 
or the maximum flux observed just before periastron 3 x 
10~ 10 erg cm -2 s _1 ). The latter choice corresponds to assuming 
that all the INTEGRAL data were taken at the X-ray maximum 
(or during the minimum assuming it is explained solely by ab- 
sorption), while the former is probably closer to the situation 
depicted by the data as the recovery had already started when 
77 Car was observed by INTEGRAL (see Fig. [TJ. 

Given the hard X-ray flux observed close to periastron by 
ISGRI, the hydrogen column density must be smaller than Nh ^ 
3 x 10 24 cm -2 , otherwise the expected flux level cannot be 
matched. Without applying this additional scaling to match the 
maximum 2-10 keV flux observed by RXTE, the extrapolation 
to the INTEGRAL range implies a lower expected flux. Less ab- 
sorption is thus needed to match the observed ISGRI flux, and 
the constraint on the hydrogen column density becomes even 
stronger, at Nh < 6 x 10 23 cirT 2 . This value is consistent with 
the highest value measured for the hydr ogen column density jus t 



Hamaguchi et al. 2007 1. 



after periastron (Nh - 4.2 x 10 23 cm -2 

The latter result also indicates that if the hydrogen column 
density during the new ISGRI observations did not increase 
above this maximum value observed during previous periastron 
passages, the flux observed by ISGRI can be explained without 
the need for an intrinsic increase in the hard X-ray emission. The 
physical conditions in the shock are expected to vary along the 
orbit. However, the average hard X-ray flux emitted seems to 
remain constant at least as close to the periastron as the ISGRI 
observations (0=3.02-3.035, corresponding to a separation of 
approximately 5-8 times the radius of the primary star). 

This lack of variability is somewhat unexpected, as the con- 
tribution of the thermal component is known to vary strongly, 
and remains fairly important (about 40%) in the first ISGRI bin. 
However, the available INTEGRAL data do not have sufficient 
statistical significance to test this. One natural reason for this 
steadiness might be that the observations were not all taken dur- 
ing the X-ray minimum (see Fig. [TJ. A possible alternative ex- 
planation of the lack of strong variability with the orbital phase 
could be that the hard X-ray emission originates, at least partly, 
from the shock between the strong winds of 77 Car and its sur- 
roundings. 

If both thermal and non-thermal components were to be af- 
fected by different absorptions, the simulations would remain 
valid provided that the hydrogen column density is larger for the 
power law than for the thermal emission. However, the upper 
limit derived would then apply to the hydrogen column density 
of the power-law component only, and that no constraints could 
be drawn for the soft X-ray part based on ISGRI observations. 

To check whether the fraction of particles injected in the 
high-energy tail is constant even at the X-ray maximum or if 
it is lower than during the rest of the orbit, future ISGRI obser- 
vations in the hard X-ray domain taken right before the X-ray 
minimum would be very useful. This would provide more con- 
straints on the evolution of the intrinsic emissivity close to pe- 
riastron. Longer observations performed in the X-ray minimum 
would also help placing a more stringent upper limit on the ab- 
sorption column density. During the next X-ray minimum (in the 
Summer 2014), 77 Car should be observable by INTEGRAL. 

On the basis of the ISGRI detection and the spectral shape, 
|Leyder et al.| ( |2008[ ) predicted the source detection in the GeV 
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energy range. The detection of 77 Car happened indeed, first with 
AGILE and subsequently with Fermi, each with comparable val- 
ues for the source flux. A graph of the spectral power flux of 
r\ Car, ranging from keV to GeV, is shown in Fig. 4 of Tavani 
et al. (2009 1. Over this wide range, the extrapolation from the 



INTEGRAL energy range up to the energy domains observed by 
AGILE and Fermi with a photon index of 1.8 (as suggested 
by these new INTEGRAL observations) suggests a flux of ap- 
proximately 5 x 10~ n ergcrrT 2 s _1 , a value compatible with the 
AGILE measurement. 

Unveiling all the mysteries related to the mechanism(s) re- 
sponsible for the X-ray minimum is likely to require observa- 
tions with the next generation of high-energy satellites, such 
as the nuclear spectroscopic telescope array (NuSTAR) or the 
Japanese project ASTRO-H. 

7. Conclusions 

A long (~90 ks) Chandra X-ray observation was analysed to 
study all sources detected in the 5-10 keV energy range around 
77 Carinae. Assuming that the INTEGRAL component were rep- 
resented by a power law (with a photon index Y of 1.8, as sug- 
gested by the new ISGRI data discussed in this paper), then this 
component would produce much more emission in the Chandra 
band than seen for any point-like source within the ISGRI error 
circle except 77 Car, provided that the hydrogen column density 
towards the ISGRI source remains modest (A^h < 6 x 10 23 cm -2 ). 
The possible contribution of the outer shell to the non-thermal 
component also remains fairly modest (below 15%). 

If the counterpart to the INTEGRAL emission were a source 
other than 77 Car, then the absorption would have to be at least of 
the order of A^h Z 10 24 crrT 2 . Such an object with a huge absorp- 
tion may be either a Compton-thick AGN serendipitously coin- 
cident with 77 Car (but the probability of this is below 0.002%), or 
a highly obscured high-mass X-ray binary (but only one source 
with the enormous hydrogen column density required is cur- 
rently known, and the hard X-ray spectral shape observed with 
ISGRI does not match that of an accreting pulsar). 

Thus, the hard X-ray emission observed with INTEGRAL is 
very likely to be associated with 77 Carinae. 

The 3<x-detection of 77 Car in the hard X-ray domain close to 
periastron can be explained by an increase in the column den- 
sity (reaching values of Nu - 6 x 10 23 cirT 2 , consistent with the 
values derived from observations during previous periastron pas- 
sages by Hamaguch i et al.|2007| l, without the need for a strong 
modification of the intrinsic emission. The situation is different 
in soft X-rays, where the emission cannot be explained only by 
the observed variations in the Nh- The energy injected in hard 
X-rays averaged over a long (about a month) timescale seems to 
be quite constant even when the companion is within the accel- 
eration region of the wind. 
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